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Recent results from short-baseline neutrino oscillation experiments and Cosmic Microwave Back- 
ground anisotropy measurements suggest the presence of additional sterile neutrinos. In this paper 
we properly combine these data sets to derive bounds on the sterile neutrino masses in the 3+1 and 
3+2 frameworks, finding a potentially good agreement between the two datasets. However, when 
galaxy clustering is included in the analysis a tension between the oscillation and cosmological data 
is clearly present. 
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I. INTRODUCTION 

In recent years, the impressive experimental discoveries 
in two fields of investigation, namely neutrino physics and 
cosmic microwave background anisotropies, have revolu- 
tionized our knowledge in particle physics and cosmology. 
Neutrino oscillations experiments have not only firmly es- 
tablished that neutrino are massive and mixed particles 
(for reviews, see e.g. Refs. [H-Hl), but have also pro- 
vided precise measurements of the three-neutrino mixing 
parameters (see the recent global fits in Refs. [1, HI). 
On the other hand, the measurements of the angular 
spectrum of the Cosmic Microwave Background (CMB) 
anisotropies (see e.g. Ref. Q) have not only fully con- 
firmed the expectations of the standard cosmological sce- 
nario but also provided a precise determination of most 
of its parameters. Moreover, with the continuous ex- 
perimental improvements, a clear interplay between neu- 
trino physics and cosmology is emerging. Neutrinos are 
indeed a fundamental energy component in modern cos- 
mology. A cosmological neutrino background is expected 
in the standard model and affects both the shape of the 
CMB and the formation of cosmological structures (see 
e.g. Ref. Q)- The recent cosmological data have pro- 
vided a clear evidence (more than 5 standard deviations) 
for the existence of the primordial neutrino background 
and have strongly constrained the absolute neutrino mass 
scale (see e.g. Ref. Q). 

However, the measurements of CMB anisotropies made 
by the ACT (Atacama Cosmology Telescope) Q and 
SPT (South Pole Telescope) experiments, when com- 
bined with the measurements of the Hubble constant 
Hq and galaxy clustering data, have provided interesting 
hints for an extra relativistic weakly interacting compo- 
nent, coined dark radiation. Parameterizing this energy 
component with the effective number of neutrino species 
iVcff, the recent data bound it to iVoff = 4.08 ±0.8 at 95% 
C.L. (see e.g. Ref. [lll - [T^ ) whereas the standard predic- 



tion for only three active neutrino species is N^s = 3.046 
[m . While this result should be taken with some grain of 
salt, since it is derived from a combination of cosmolog- 
ical data and some tension does exist between the data 
(see e.g. Ref. [l^) it is anyway interesting since a fourth, 
or fifth, neutrino species seems also suggested by short- 
baseline (SBL) oscillation experiments. The appearance 
and disappearance data of several SBL experiments can 
be explained by the mixing of the three active neutri- 
nos with one or two additional sterile neutrinos in the 
so-called 3+1 and 3+2 models (see Refs. [l3-il). 

This work is aimed to determine the masses of the 
sterile neutrinos in 3+1 and 3+2 models using data from 
SBL experiments and recent cosmological data and check 
if the results are mutually compatible. Finally, we com- 
bine the bounds from the two different analyses to have a 
joint probability for the masses of sterile neutrinos. The 
paper is organized as follows: in Sec. and in Sec. IIIII 
we present the data sets we make use of, the method we 
adopt to analyze them and the results we obtain regard- 
ing the SBL experiments and in the cosmological context, 
respectively; in Sec. llVl thc joint analysis method and re- 
sults arc shown; finally we summarize our conclusions in 
Sec.lVl 



II. NEUTRINO OSCILLATIONS ANALYSIS 

The short-baseline neutrino oscillation analysis is per- 
formed following Refs. [H-i^. 

We consider 3+1 and 3+2 neutrino spectra in which 
i/g, Vfj,, Vt are mainly mixed with v\, vi, v^, whose masses 
are much smaller than 1 eV and there are one or two ad- 
ditional massive neutrinos, 1^4 and 1^5, which are mainly 
sterile and have masses of the order of 1 eV. Short- 
baseline oscillations are generated by the large squared- 
mass differences Amli and Arn^i, with: 



(1) 
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The small squared-mass differences Arn^i and Am'^i 
which generate, respectively, solar and atmospheric neu- 
trino oscillations (see Refs.[l|-[3]) have negligible effects 
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3+1 


3+2 


Xmin 


142.1 


134.1 


NDF 


138 


134 


GoF 


39% 


48% 


Am|i [eV^] 

|f/e4|' 


1.62 


0.89 


0.035 


0.018 


Amii [eV^] 


0.0086 


0.015 




1.61 




0.022 






0.0047 






1.577r 
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TABLE I: Values of Xmim number of degrees of freedom 
(NDF), goodness-of-fit (GoF) and best-fit values of the mix- 
ing parameters obtained in our 3+1 and 3+2 fits of short- 
baseline oscillation data. 



in SBL oscillations and are ignored in the following. 
The two heavy neutrino masses 7714 and m.g which are 
probed by cosmological data are simply connected to the 
squared-mass differences relevant for SBL oscillations by: 



1714 



7715 



'51 



(2) 



We fit the data set of short-baseUne neutrino oscillation 
experiments corresponding to the GLO-HIG analysis in 
Ref. [2^, in which the low-energy MiniBooNE neutrino 
[23} and antineutrino [23 - l26j data corresponding to the 
so-called "MiniBooNE low-energy anomaly" are not con- 
sidered, since they induce a strong tension between ap- 
pearance and disappearance data (see the discussions in 
Refs. [il,!!^]). We made the following two improvements 
with respect to the analysis presented in Ref. (20l |: 

1 . We used the reactor neutrino fluxes presented in the 
recent White Paper on light sterile neutrinos j27j . 
which update Refs. [2^ IsO]- The new fluxes are 
about 1.3% larger than those we used before, which 
were taken from the reactor antineutrino anomaly 
publication (29| . 

2. We replaced the KamLAND bound on |(7e4|^ with 
a more powerful constraint obtained from solar 
neutrino data f3ll - l33j . Taking into account the 
recent measurement of |J7e3p in the Daya Bay 
[3^ and RENO [s^ reactor neutrino experiments 

|{7e3f_= sin^iJia = 0.025 ± 0.004), from Fig. 1 of 



Ref. [33| we inferred the approximate upper bound 
|J7e4p = sin^du < 0.02 at la (see Ref. [36|]). 

In our analysis of SBL neutrino oscillation data we ap- 
ply first the standard method. The minimum value 
of x^i the number of degrees of freedom, the goodness- 
of-fit and the corresponding best-fit values of the oscil- 
lation parameters are presented in Tab. HI The results 
concerning the 3+1 and 3+2 fits are similar to those re- 
ported, respectively, in Ref. [l^l for the GLO-HIG case 
and Ref. |18j , with small variations due to the considera- 
tion of different data sets. From Tab. U we can see that in 
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FIG. 1: 



Allowed regions in the sin^ 2i?e/j-A7Ti4i, sin^ 2i9e 



A?7i4i and sin 2i9p^-A7n4i planes obtained from the global fit 
of short-baseline neutrino oscillation data in the 3+1 scheme 



using the standard x method, 
by a cross (see Table. |l}. 



The best-fit point is indicated 



both the 3+1 and 3+2 frameworks the global goodness- 
of-fit is satisfactory. 

The allowed regions of Amli versus the effective SBL 



oscillation amplitudes sin 21}^^ , 

9a0 = 4|[/Q4p|[/^4 



sin^ 2dpp and sin^ 2?9,, 



(with sin^ 279ct,9 = 4|[/q4P|[/^4|^) are shown in Fig. [T] 



(-) 



These regions are relevant, respectively, for ly^j, 

— >Ve andV^ -^'i^li oscillation experiments. They are 
more similar to those shown in Fig. 3 of Ref. [l^l than 
the region presented in Ref. [11] , because the larger reac- 
tor antineutrino fluxes used in this analysis increase the 
reactor antineutrino anomaly, leading to a larger value 
of |t/e4p, which tends to cancel the effect of the solar 
neutrino constraint. 

The allowed regions in the Am|]^-A77igj^ plane obtained 
in the 3+2 analysis are shown in Fig. [2] One can see 
that the allowed regions are similar to those presented 
in Fig. 9 of Ref. [l^l: with small variations due to the 
different considered data sets. 

Since we want to perform a combined analysis of SBL 
oscillation data and cosmological data and the cosmo- 
logical analysis is performed with the Bayesian method, 
we have also analyzed the SBL oscillation data with a 
Bayesian approach. We assumed the sampling distribu- 
tion of the data D: 



(3) 



where M is the model (M = 3 + 1 or A/ = 3 + 2), 6* a/ is 
the corresponding set of oscillation parameters (listed in 
Tab.|T| and x^{D,9m) is the corresponding function. 
The sampling probability is called "likelihood" when con- 
sidered as a function of the parameters of the model. In 
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3+2 

68.27% C.L. (1o) 
90.00% C.L. 
95.45% C.L. (2o) 
99.00% C.L. 
99.73% C.L. (3c) 




FIG. 2: Allowed regions in the Amlj-Amli plane and cor- 
responding marginal Ax^'s obtained from the global fit of 
short-baseline neutrino oscillation data in 3+2 schemes using 
the standard method. The best-fit point is indicated by a 
cross (see Table. [l|. 



each of the two models, we calculated the posterior prob- 
ability distribution of the oscillation parameters using 
Bayes' theorem: 



pieM\D,M) 



p{D\eM,M)p{9M\M) 
p{D\M) 



(4) 



where p{D\M) is easily calculated as a normalization 
constant. We assumed a flat prior distribution in the 
logarithmic space of the oscillation parameters, except 
for the CP-violating phase rj in the 3-1-2 spectrum (see 
Rcf. [l^) for which we used a linear scale in the inter- 
val [0,27r]. For log(ATO|i/eV^) and log(Am|i/eV^) we 
considered the range [—1,1]- For log|C/e4p, log|t/^4p, 
log |t/e5p, log |?7;i5p we considered the range [—4,0]. 

Since we are interested in combining the results of the 
analysis of SBL oscillation data with that of the cosmo- 
logical data, where the only shared parameters are the 
neutrino masses in Eq.([2]), we calculated the marginal 
posterior probability distributions of the squared-mass 
differences by integrating the posterior probability distri- 
bution over the other oscillation parameters taking into 
account the scale of the flat prior. For example, in the 
3+1 model: 



p(logAm|i|Z?,3 + l) = 
X p(log(Am^i), log |[/e4p, log \U^4\^\D, 3 



d\0g\Uei\^ d\og\U^,i\^ 



(5) 



In this way, we obtained the posterior probability distri- 
bution of Amli in the 3+1 spectrum plotted in Fig. [3] 
(thick green line exhibiting several sharp peaks) and the 
allowed regions in the Amlj^-Arngj^ of the 3+2 spectrum 



shown in Fig. |4l Comparing with Fig. [2l one can see 
that the Bayesian allowed regions are wider than those 
obtained with the method. The difference is due to 
the different method of marginalization with respect to 
the other mixing parameters (mixing angles and CP- 
violating phase): in the method one considers only 
the minimum of the in the range of each marginal- 
ized parameter, whereas in the Bayesian method one 
must integrate the posterior probability density over the 
marginalized parameter space. Since the data do not con- 
strain much the values of the marginalized parameters 
(see Figs. 10-12 of Ref. [13 )i the Bayesian integration 
gives significantly different results from the x^ marginal- 
ization. The allowed vertical bands with constant value 
of Amli are due to the fact that one can have a compa- 
rable fit for any value of Am^^^ and negligible \Ue5\ and 
C/^isI, which is effectively equivalent to a 3+1 framework. 
The same applies to the allowed horizontal bands with 
constant value of Amli . 



III. COSMOLOGICAL ANALYSIS 

The cosmological analysis is performed in two different 
steps: first by analyzing CMB-only data and then by fur- 
ther adding data from large scale structure and priors on 
the Hubble parameter. The CMB analysis is performed 
by employingthe following datasets: WMAP7 Q, ACT 
Q and SPT [1^ . The large scale structure analysis makes 
use of information on dark matter clustering from the 
matter power spectrum extracted from the SDSS-DR7 
luminous red galaxy sample |37| . Finally, the Hubble pa- 
rameter prior we use is based on the latest Hubble Space 
Telescope observations [38l |. 

Wc analyze datasets up to €max = 3000. The analy- 
sis method we adopt is based on the publicly available 
Markov Chain Monte Carlo (MCMC) package CosmoMC 
{39} with a convergence diagnostic done through the Gel- 
man and Rubin statistic. 

We sample the following six-dimensional standard set 
of cosmological parameters, adopting flat priors on them: 
the baryon and cold dark matter densities 17b /i^ and 
flch'^, the ratio of the sound horizon to the angular di- 
ameter distance at decoupling 9, the optical depth to 
reionization r, the scalar spectral index ns and the over- 
all normalization of the spectrum As- We account for 
foregrounds contributions including three extra ampli- 
tudes: the SZ amplitude, the amplitude of clustered point 
sources, and the amplitude of Poisson distributed point- 
sources. Wc consider purely adiabatic initial conditions 
and we impose spatial flatness. In this work both active 
and sterile neutrinos are assumed to be fully thermalized 
(for the non thermal case see, e.g., Ref. j4(j|). 

The aim of this paper is to specifically test 3+1 and 
3+2 neutrino mass models, by means of a joint analysis of 
both cosmological and SBL experiments data. Therefore, 
contrary to the typical approach (see e.g. Ref. (4l| - |43j ). 
in the cosmological analysis we do not let the effective 
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number of relativistic degrees of freedom A^off to vary as a 
free parameter, instead we fix it at the values iVoff = 3 + 1 
or iVefi = 3+2 for the 3+1 and 3+2 schemes, respectively. 
This is consistent with the assumptions done in the oscil- 
lation analysis and with the hypothesis of cosmological 
full thcrmalization of all neutrino states (includin g th e 
sterile ones; see the recent discussions in Rcfs. (i^. |45|). 
Consistently to the analysis of Section [TTl we fix the three 
active neutrinos to be massless and we allow the sterile 
neutrinos to have masses which vary as additional free 
parameters. Since we are interested to sample the joint 
sensitivity of cosmological and SBL neutrino data on the 
sterile-neutrinos mass parameters, in the cosmological 
analysis we do not employ the neutrino mass fraction 
fu (as it is usually done), but instead we sample directly 
log Anili and log Ato|]^ . This implies a fiat prior on those 
parameters. 

Before attempting a joint analysis with the SBL data, 
which have been presented in the previous Section, we re- 
port in Tab.mithc constraints on the cosmological param- 
eters using CMB-only data and CMB data plus SDSS in- 
formation together with the HST prior, and assuming: a 
3+1 model with three massless active neutrinos and one 
massive sterile neutrino; a 3+2 model with three mass- 
less active neutrinos plus two massive sterile neutrinos. 
The 95% C.L. mass bounds on the sterile neutrinos is 
2.88 eV for the 3+1 scheme, while for the 3+2 model the 
bound on the sum of the masses of the two additional 
sterile neutrinos is 2.48 eV, both of them share a 2a up- 
per limit of about 1.24 eV, when CMB-only data are 
used. These bounds drastically improve when also SDSS 
data and the HST prior are included in the analysis (see 
Ref. 53), reaching the value of 0.73 eV for the 3+1 case 
and about 1 eV for the 3+2 case. Both the 3+1 and 3+2 
schemes are statistically well acceptable, with no notice- 
able preference in the minimal x^- The only visible (and 
expected) difference between the 3+1 and 3+2 schemes 
is that 2 additional neutrinos require a larger value of the 
dark matter abundance ilch^, to compensate a delay of 
the equivalence time, which would instead be induced by 
the presence of an additional light degree of freedom in 
the 3+2 case [131 . The correction due to non degener- 
acy between the mass of the first and the second sterile 
neutrino in the 3+2 model is of the order of precision of 
present numerical codes and so undetectable using only 
the present cosmological data (CMB and matter power 
spectra). Moreover the degeneracies with other cosmo- 
logical parameters makes the detection of the neutrino 
mass differences impossible at the state of art (see Ref. 

M)- 

Fig. |3] shows the marginal posterior probability of the 
cosmological Bayesian analysis for the 3+1 case, com- 
pared with the results of the SBL study. The blue line 
exhibiting a broad peak stands for the analysis of the cos- 
mological data alone and the left panel refers to CMB- 
only data, while the right panel refers to the CMB data 
implemented with SDSS and HST information. The two 
panels of the figure show how the inclusion of SDSS and 



HST information is relevant to set the more stringent 
constraint on the cosmological upper bound on the neu- 
trino mass. The shaded regions refer to the 95% C.L. 
coverage of the probability distribution, from which the 
bounds on TO4 of Table [IT] are derived. When com- 
pared with the SBL analysis and its 95% C.L. mass in- 
tervals (three slightly discontinued ranges in the inter- 
val 0.93 eV < m4 < 1.45 eV and a higher mass range 
2.29 eV < TO4 < 2.59 cV), with a best fit at TO4 = 1.27 eV, 
we notice that CMB-only and SBL oscillation data are 
well compatible among them, with a significant overlap 
of the corresponding 95% C.L. regions. The 95% C.L. 
cosmological upper bound TO4 < 2.88 eV disfavors the 
higher mass SBL solution, while is perfectly compatible 
with the lower SBL mass ranges. The combination of the 
cosmological and SBL datasets will therefore produce a 
clean allowed interval, as shown in the next Section. In- 
stead, when SDSS and HST information are included in 
the analysis, SBL oscillations and cosmological data are 
in tension, with no overlapping 95% C.L. 

The analysis for the 3+2 scheme is shown in Fig. UJ 
where C.L. regions in the Amli-Arn^i plane are re- 
ported. The SBL allowed regions clearly show a pref- 
erence for at least a non-zero neutrino mass (7715 with 
our choice of hierarchy in neutrino masses) and a global 
preference for 7714 = 0.95 eV and 771,5 = 1.27 eV. The 
cosmological data instead provide upper limits on both 
sterile neutrino masses, with no clear preference for non- 
zero values. CMB-only data (left panel) are well compat- 
ible with SBL results, with the 95% C.L. upper bound of 
the cosmological analysis consistent with the correspond- 
ing 95% C.L. regions of the SBL analysis and its global 
best-fit point (7774 = 0.95 eV and = 1.27 eV). Also 
in the 3+2 case, the inclusion of SDSS and HST data 
produces tension between SBL and cosmological analy- 
ses, as is manifest in the right panel of Fig. 21 where 
only a partial overlap at the 3cr C.L. is present. Fig. |4] 
clearly shows that the whole set of cosmological data will 
be instrumental in significantly reducing the degeneracy 
of the allowed solutions of the SBL analysis when the 
joint analysis will be attempted in the next Section. 



IV. COMBINED ANALYSIS 

The combined analysis of the SBL oscillation data 
and the cosmological observations has been performed 
by merging the corresponding posterior probabilities. 
Since the only relevant parameters common to both sec- 
tors are the sterile neutrino masses 7774 ~ y^Am^ and 
7715 ~ Am|j^ we can define a marginal posterior proba- 
bility for the joint analysis by directly multiplying the 
SBL and cosmological marginal posterior probabilities 
relative to the parameter of interest. For example, in 
the 3+1 case, denoting by Dq and Ds the cosmological 
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3+1 CMB only 


3+2 CMB only 


3+1 CMB+SDSS+HST 


3+2 CMB+SDSS+HST 


T 

Ho 

ris 


0.0224 ± 0.0004 
T^5 + 007 
0.085 ±0.014 

71.5 ±3.6 
0.970 ±0.015 
3.21 ±0.05 


0.0226 ± 0.0004 
1 56 + 009 
0.087 ±0.015 

73.6 ±4.4 
0.985 ±0.016 
3.20 ±0.05 


0.0224 ± 0.0004 
133 + 004 
0.084 ± 0.014 

73.1 ± 1.6 
0.977 ±0.010 
3.19 ±0.04 


0.0226 ± 0.0004 
1 56 + 004 
0.086 ±0.014 

74.6 ±2.0 
0.990 ±0.010 
3.19 ±0.04 


Em (eV) 


< 2.88 


< 2.48 


< 0.73 


n (+0.45) 
"•■^°-0.13 (-0.42) 


■ 

'^niin 


7529.5 


7532.2 


7578.5 


7581.1 



TABLE II: MCMC estimation of the cosmological parameters from the analysis of CMB-only data and from CMB data plus 
matter power spectrum information (SDSS) and a prior on Ho (HST), in the case of three massless active neutrinos and one 
massive sterile neutrino (3+1 scheme) and assuming 3 massless active neutrinos plus 2 massive sterile neutrinos (3+2 scheme). 
Neutrino mass upper bounds are reported at the 95% C.L., unless for the 3+2 CMB+SDSS+HST case where we quote the 
best-fit value together with the 68% (95%) C.L. interval. 




FIG. 3: Marginal posterior probabilities obtained with a Bayesian analysis for Amli in the 3+1 scheme. The thick [green] 
solid line exhibiting several sharp peaks (the same in the two panels) refers to the analysis of the short-baseline oscillation data 
alone. The blue line exhibiting a broad peak stands for the analysis of the cosmological data alone: CMB-only data for the left 
panel, CMB data implemented with SDSS and HST information for the right panel. In all cases, the shaded regions refer to 
the 95% coverage of the probability distribution. 



and SBL data we have^: 

p{\ogAml,\Dc+s,5 + l) (x (6) 
p{\ogAml^\Dc,S + l) X p{\ogAml^\Ds,3+l), 

where the SBL probabihty is the one defined in Eq. ([5]) 
and the cosmological probability is the one used in the 



Am|j in both the SBL 



^ Since we assumed a flat prior for 9 = 
and cosmological analyses, using Bayes' theorem ijJJ we have 

p{e\Dc+s)'xp{Dc+s\e) = p{Dc\e)p{Ds\e)^p{e\Dc)p{e\Ds) 



analysis of the previous section and obtained through 
CosmoMC. 

The combined analysis for the 3+1 scheme is shown 
in Fig. [51 As usual, the left panels refers to the case 
of CMB-only data in the cosmological sector, while the 
right panel adds SDSS and HST datasets. The horizon- 
tal dashed lines identify the credible intervals at 68.27%, 
90.00%, 95.45%, 99.00% and 99.73% C.L. In the case 
of CMB-only data, the inclusion of the cosmological in- 
formation to the SBL analysis disfavors the higher mass 
SBL solution around 2.4 cV but maintains the lower mass 
95% C.L. allowed intervals (0.90 cV < < 1.46 cV) 
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FIG. 4: Allowed regions in the Amii-Amli plane obtained with a Bayesian analysis. The "boxy" regions (the same in the 
two panels) refer to the global analysis of the short-baseline oscillation data and are relative to the following confidence levels 
(from the innermost to the outermost region): 68.27% (red), 90.00% (light blue), 95.45% (green), 99.00% (brown) and 99.73% 
(dark blue). The arc-shaped solid lines refer to the analysis of the cosmological data: the left panel stands for the CMB-only 
dataset, while the right panel refers to the inclusion of the SDSS information and HST prior to the CMB data. The different 
lines refer to the following confidence levels (from the lower curve to the upper curve, in each panel): 68.27% , 90.00%, 95.45%, 
99.00% and 99.73% . 
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FIG. 5: Marginal posterior probabilities obtained with a Bayesian analysis for Amli in the 3+1 scheme, for the joint analysis 
of cosmological and short-baseline data. Left panel: the cosmological analysis employs CMB-only data. Right panel: the 
cosmological analysis adds SDSS and HST information to the CMB data. The horizontal dashed lines identify (from the lower 
curve to the upper curve, in each panel) the credible intervals at 68.27%, 90.00%, 95.45%, 99.00% and 99.73% C.L. 



and (2.27 cV < 7714 < 2.51 eV) and the best-fit solution 
(to4 1.27 eV). When SDSS and HST information is 
added to the analysis, the allowed interval of the global 
analysis shifts down to lower values of the sterile neu- 
trino mass, due to the more stringent bound from the 
cosmological sector. The 95% C.L. mass range becomes 
0.85 eV < 7714 < 1.18 eV, and the best fit shifts down to 



TO4 = 0.93 eV. 

The combined analysis for the 3+2 scheme is shown in 
Fig. ini again for the case of CMB-only data (left panel) 
and for the further inclusion of SDSS and HST data (right 
panel). The global results are that at least one sterile 
neutrino needs to be massive, with a mass of the order of 
1 cV (tos with our choice of hierarchy) , while the second 
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FIG. 6: Allowed regions in the Amli-Amii plane obtained with a Bayesian approach, for the joint analysis of short-baseline 
and cosmological data. The different regions (as in Fig. |4]) refer to the following confidence levels (from the innermost to the 
outermost region): 68.27% (red), 90.00% (light blue), 95.45% (green), 99.00% (brown) and 99.73% (dark blue). Left panel: 
SBL data plus the CMB-only dataset. Right panel: SBL data plus CMB, SDSS and HST data; in this case only 68.27% (red), 
95.45% (green) and 99.73% (dark blue) C.L. are reported. 



sterile neutrino can be massless. The marginalized 95% 
intervals for the two neutrino masses are: < 2.51 eV 
and 0.86 eV < < 3.16 eV when CMB-only data are 
considered; 7714 < 0.70 cV and 0.67 cV < < 1.35 cV 
for the full analysis which includes also SDSS and HST. 



V. CONCLUSIONS 

Measuring the number and the mass of sterile neutri- 
nos is one of the most interesting challenges both in cos- 
mology and in neutrino physics. The existing cosmologi- 
cal data indicate that the energy density of the Universe 
may contain dark radiation composed of one or two ster- 
ile neutrinos, which may correspond to those in 3-1-1 or 
3-1-2 models which have been invoked for the explanation 
of short-baseline neutrino oscillation anomalies. We have 
performed analyses of the cosmological and SBL data in 
the frameworks of both the 3-1-1 and 3-1-2 models. Then 
we have compared the results obtained with the same 
Bayesian method, to figure out if the indications of cos- 
mological and SBL data are compatible. 

At the state of art, cosmological data are sensitive to 
the sum of neutrino masses, for which they give an up- 
per limit at the scale of about 1 cV. Hence they do not 
allow us to resolve the degeneracy between the mass of 
the first and the second sterile neutrino in a 3-1-2 model, 
although in the numerical calculation we leave them as in- 
dependent parameters. Instead, short-baseline neutrino 
oscillations have a completely different parameterization 
and in the 3-1-2 model the degeneracy between the two 
square mass differences Am|j^ and Ato|]^ is broken. 



The results of our analysis show that the cosmological 
and SBL data give compatible results when the cosmo- 
logical analysis takes into account only CMB data. But 
if the information on the matter power spectrum com- 
ing from galaxies surveys are also considered there is a 
tension between the sterile neutrino masses needed to 
have SBL neutrino oscillations and the cosmological up- 
per limit on the sum of the masses. 

The combined analysis of cosmological and SBL data 
gives an allowed region for in the 3-1-1 scheme around 
1 eV. In the 3-f2 scheme, the cosmological data reduce 
the allowance of the second massive sterile neutrino given 
by SBL data, leading to a combined fit which prefers the 
case of only one massive sterile neutrino at the scale of 
about 1 eV. 

In conclusion, our analysis shows that cosmological 
data are marginally compatible with the existence of one 
massive sterile neutrino with a mass of about 1 eV, which 
can explain the anomalies observed in SBL neutrino oscil- 
lation experiments. The case of massive sterile neutrinos 
is less tolerated by cosmological data and in any case the 
second sterile neutrino must have a mass smaller than 
about 0.6 eV. 
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